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We develop a systematic field theoretic description for the roughness correction to the Casimir free 
energy of parallel plates. Roughness is modeled by specifying a generating functional for correlation 
functions of the height profile, the two-point correlation function being characterized by the variance, 
a 2 , and correlation length, I, of the profile. We obtain the partition function of a massless scalar 
quantum field interacting with the height profile of the surface via a 5-function potential. The 
partition function of this model also is given by a holographic reduction to three coupled scalar 
fields on a two-dimensional plane. The original three-dimensional space with a flat parallel plate at 
a distance a from the rough plate is encoded in the non-local propagators of the surface fields on its 
boundary. Feynman rules for this equivalent 2 + 1-dimensional model are derived and its counter 
terms constructed. The two-loop contribution to the free energy of this model gives the leading 
roughness correction. The absolute separation, a e ff, to a rough plate is measured to an equivalent 
plane that is displaced a distance p oc a 2 /£ from the mean of its profile. This definition of the 
separation eliminates corrections to the free energy of order l/a* ff and results in a unitary model. 
We derive an effective low-energy theory in the limit < < a. It gives the scattering matrix and 
equivalent planar surface of a very rough plate in terms of the single length scale p. The Casimir 
force on a rough plate is found to always weaken with decreasing correlation length I. The two-loop 
approximation to the free energy interpolates between the free energy of the effective low-energy 
model and that of the proximity force approximation - the force on a very rough Dirichlet plate 
with a > 0.5£ being weaker than on a flat plate at any separation. 
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I. INTRODUCTION 

Casimir orig inally[l] obtained the force due to electromagnetic zero-point fluctuations between two large ideal 
parallel metallic flat surfaces at vanishing temperature. His approach was soon generalized to dielectric surfaces^, 
finite temper atureQ, and experimentally more accessible geometries Q. The influence of surface roughness was 
considered only much later[^, perhaps because this correction was insignificant in early Casimir experiments. 
When Casimir forces were accurately measured with atomic force microscope techniques [7| at plate separations of 
only a few hundred nanometers, corrections caused by the roughness of the plates could no longer be ignored. They 
are even more important at the small separations and higher accuracy of more recent experiments^^. An increasing 
amount of experimental Q and theoretical [9j-[l5j effort has since been devoted to understanding roughness effects. 
Currently the only ri gor ous non-perturbative approach is the proximity force approximation(PFA) (and some recent 
modifications thereof [15| ) ■ It is accurate when the correlation length £ of the profile greatly exceeds the separation 
a of the plates [13, EH- Most other approaches consider perturbative corrections to the Green's function in powers of 
a/a. The limit of extremely rough plates with a»i was first obtained using methods of stochastic calculus |5;j. 

For stochastic roughness all perturbative calculations @, H, E2 show an increase in magnitude of the Casimir energy 
and force with decreasing correlation length £, approaching the PFA for £ ^> q[I3|. For a massless scalar field this 
trend is shown by the dashed curves in fig. [5j It is qualitatively similar for the electromagnetic case [l, EH. From the 
point of view of the multiple scattering expansion to the Casimir energy this strengthening of the force is not intuitive 
and in fact is physically untenable in the limit a><: the back-scattering, i.e. echo, from a rough plate generally 
decreases with increasing roughness. The reflection coefficient for scattering off a rough surface therefore ought to be 
reduced in magnitude from that for perfect reflection off a Dirichlet (or ideal metallic) plate. Since the Casimir energy 
is related to the trace of the Green's function one semi-classically expects a reduction in magnitude of the Casimir 
energy (and force) for a rough plate and not an enhancement. 

This qualitative argument becomes rigorous upon using the GTGT-formula 16J for the Casimir energy due to a 
massless scalar field for a stochastically rough but otherwise ideal (Dirichlet) plate and a perfectly smooth parallel 
flat plate. When the plate separation a is large compared to the correlation length I of the profile one recovers 
translational invariance and the reflection coefficient becomes diagonal in the transverse momentum. The Casimir 
energy per unit area, £(a), in this case is given by a dimensionally reduced atat- formulaflil [l7j . 

r°° K 2 dK 

£{a)= J "2^ ln(1 ~ W(*)e~ 2aK ) , (1) 

where the reduced reflection matrices, trough and i, for back-scattering off the two parallel plates for a massless scalar 
field are functions of the wave-number k only. The reduced scattering matrix t for the flat Dirichlet plate is t = tr> = I 
and unitarity demands that |t r ou g h(K)| < 1 = to for reflection off the rough plate. Inspection of Eq.(H|) then implies 
that £ (a) for interaction with a stochastically rough Dirichlet plate should be reduced compared to the Casimir energy 
per unit area for interaction with a flat Dirichlet plate. This argument does not hold in the regime of large correlation 
length of the PFA and is rigorous only in the limit £ — > in which translational invariance is recovered. A problem 
arises because the trend of the perturbative analysis is the opposite and gives a Casimir energy that increases in 
magnitude beyond all bounds for £ — > 0. 

The perturbative analysis for electromagnetic fields of[j| H, EH predicts large increases that appear not to be 
supported by experiment 0, Q . Some experiments @ with gold coatings described by I ~ 35nm, a ~ 5nm clearly are 
in the rough regime with a ^> £ for separations a > fOOnm. 

This investigation was partly motivated by a desire to reconcile the strengthening of the force with increased 
roughness observed in all perturbative calculations d, H, H, EH, EH with the weakening demanded by unitarity and the 
multiple-scattering formalism. [We will see that the two approaches differ in the definition of the plate separation.] 
However, the following field theoretic description goes beyond the original objective of a more rigorous and non- 
perturbative description of roughness effects. It provides a framework for a consistent loop expansion and includes 
temperature effects. The loop expansion is uniform in both small parameters a /a and a/£ and the field theory is 
interesting in itself. It is holographic [1 8] in the sense of being equivalent to a lower-dimensional field theory on the 
two-dimensional plane that is a boundary of the original space. The existence of a distant plate in an extra dimension 
is encoded by non-local propagators of the surface fields in the latter model. Roughness corrections to Casimir energies 
in this sense can be described by a table-top brane theory. 



II. THE GENERATING FUNCTIONAL OF ROUGHNESS CORRELATIONS 

We consider the standard Casimir configuration of two parallel flat plates at an average separation that is much 
less than their transverse dimensions [if. A Cartesian coordinate system with z-axis normal to the plates is used to 
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(3) 



describe this system. The profile function h(x) associated with a plate at the mean height (z) = a gives the precise 
position of this surface as a function of the (two) transverse coordinates!]] x = (x, y), 

z(x) = a + h(x). (2) 

We assume the profile of the plate is without enclosures and that h(x) is a single- valued function. It nevertheless 
is often more practical to characterize a rough plate by just a few low-order correlation functions of its profile than 
by the profile itself. For sufficiently large plates, a description in terms of (all) correlation functions in fact is exact. 
The formalism developed below may, in principle, also be applied to plates whose profile is known precisely. 

The n-point correlation functions of the profile h(x) for a plate of (large) area A are the averages, 

D 1 = {h(xi)) :=A -1 / h(x + xi)dx 

J A 

L> 2 ( Xl -x 2 ) = (/i(xi)/i(x 2 )) := A^ 1 I h(x + x 1 )h(x + x 2 )dx 

J A 

.D n (xi - x 2 , . . . ,x n _i - x n ) = (/i(xi) . . . h(x n )) := A^ 1 / h(x + xi) . .. h(x + x n )dx . 

J A 

We have here assumed that the plate is large enough for its boundary to be ignored and have used translational 
invariance to assert that the correlation functions in this case depend only on differences in the transverse coordinate^. 
Isotropy of the profile yields further restrictions; the 2-point correlation function D 2 in this case depends only on the 
distance between the two points. We assume that the profile and therefore all n-point correlation functions of Eq.([3|) 
can, at least in principle, be measured when the plate is far removed from any other object. The mean position 
(z) = a of the plate in Eq.Q is fixed by requiring that 

£>i = (h(x)) = 0. (4) 

It is convenient to collect all correlation functions of Eq. @ in a single generating functional Zh [a] , 

Z h [a} = '^2—^ if a(x 1 )a(x 2 )...a(x n )D n (x 1 ,...,x n )dx 1 dx 2 ...dx n , (5) 

and directly model Zh[a] instead of individual correlation functions. With the restriction of Eq.Q, the simplest 
model for a rough plate is entirely determined by the 2-point correlation D 2 of the profile. The generating functional 
for such a (quadratic) Gaussian model is of the form, 

4 2) H = exp±{a| J D 2 |a} , (6) 

with, 

{a\D 2 \a} := jj a(x)D 2 (x - y)a(y)dydx . (7) 

In general, Eq.(j6|) just gives the leading term in a cumulant expansion of Zh ■ Stochastic roughness is fully described 
by the covariance of the profile and a Gaussian model by definition is exact in this case. A Gaussian model will also 
suffice to extract corrections to the free energy to leading order in the roughness profile. To leading order in the 
variance a 2 even the effect of corrugated profiles /i^(x) = crsin(a;j;) can be described by a Gaussian model. But the 
four-point correlation in this case is only half of what the Gaussian model predicts, 

2 

D%{x-y) = y cos(lu{x - y)) but , (8) 

^(x 1 ,x 2 ,x 3 ,x 4 ) = ±{D%(x 1 -x 2 )D%{x 3 -x 4 )+ 

+ D%{x x - x 3 )D%(x 2 - x 4 ) + D£(xi - x 4 )D%(x 2 - x 3 )) . 



1 The normal direction is distinguished and bold-faced letters describe two-dimensional transverse vectors in the following 

2 For exact translational invariance, the finite flat plate should be replaced by a two-dimensional torus of area A. 
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To correctly describe effects due to a periodic profile to order a 4 thus requires inclusion of a 4 th order cumulant. Note 
that correlations of periodic profiles are not positive definite and have no probabilistic interpretation. 

The mathematical basis for a field theoretic approach to roughness is that any analytic functional F[h] of the profile 
h(x) with translation invariant coefficients can be evaluated using the generating functional Zy\a\. To see this, first 
consider the evaluation of a monomial in the Taylor expansion of F[h] for small profiles, 



dx x dx 2 . . . dx n F n (x 1 -x 2 ,...,X„_i - x n )h(x 1 )h(x 2 ) . . . h(x n ) = 

i- J dx jj dx x dx 2 . . . dx n F n (xi - x 2 , . . . ,x n _i - x n )ft(x + xi)ft(x + x 2 ) . . . h(x + x n ) 



A 

= jj dx x dx 2 . . . <2x„F„(xi - x 2 , . . . , x„_i - x„)D„(xi - x 2 , . . . , x n _i - x„) 

6 6 6 



(9) 



c?xi<ix 2 . . . dx n F n (xi - x 2 , . . . , x n _i - x, 



<5a(xi) 6a(x 2 ) 6a(x n ) 



Zh[a] 



a=0 



The first equality in Eq.© is due to the translational invariance of the coefficient functions F n [but assumes no 
regularity of the profile h(x) itself]. No further assumptions are required and Eq.© holds for any profile on a 
sufficiently large plate. The third line in Eq.© implies that the result is proportional to the area A. Assuming that 
all coefficient functions F n in the Taylor expansion of F[h] are translation invariant and that the expansion converges 
for the particular profile, Eq.© implies that one may formally evaluate F[h] by the functional derivative, 



F[h] = F[—]Z h [a] 
da 



a=0 



(10) 



It remains to (approximately) determine the dependence of the partition function of a massless scalar field on the 
profiles of parallel plates. 



III. DEPENDENCE OF THE FREE ENERGY OF A MASSLESS SCALAR ON THE PROFILES OF 

PARALLEL PLATES 

This problem has to some extent been addressed in the calculation of lateral Casimir forces for corrugated plates (l7l 
[l9j . but regular, one dimensional corrugations are very special. Lateral Casimir forces are automatically finite and 
vanish if one of the two parallel plates is flat. Here we are interested in the influence of profiles on the normal Casimir 
force. The physical interpretation and consistent subtraction of (divergent) contributions due to the curvature of the 
profile function is one of our main objectives. 

As in ref.[i3], we model the interaction of the n th thermal mode[2l| of a massless scalar quantum field <p n (x,z) 
[corresponding to Matsubara frequency = 2imT] with two semitransparent parallel plates by (5-function potential^, 

Vint(x, z) = \6{z - h(x) -a) + X6(z - h(x)) . (11) 

Here h, h <C a/2 are the profiles of the two surfaces at average positions (z) = a and (z) = 0, respectively. A and 
A > are the corresponding coupling constants of canonical length dimension a -1 . The limit A — > oo (or A — > oo) 
suppresses tunneling through, and enforces Dirichlet boundary conditions on, the corresponding surface. For finite 
coupling the plate is semitransparent. 

Although this scalar model appears far removed from reality, it is sufficiently simple to analyze thoroughly and does 
exhibits some features encountered in the clcctrodynamic case. It in particular essentially describes the thin plate 
limit of the electric contribution to the Casimir force (22J- Since we are not primarily interested in lateral forces [l7l[T9|. 
we consider the case where only one of the plates is rough by setting h — 0. This restriction does not qualitatively 
affect the approach but greatly simplifies the model. Expanding the interaction of Eq. (|ll[) for h(x) a, one arrives 



3 We use natural units h = c = ks = 1 throughout. 
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at the interaction Hamiltonian, 



-Hint [h, < 



/ rfxi[A^(x,a + / l (x)) + A^(x,0)] 

n J 

[h] + J2 ffffi fen] + H$ [h, 4> n } + [h, d> n ] - 



with 



(0), 



H; 



= d. 



~4> 2 (x,a) + ~^(x,0) 



A 



H^[h] = 



h m (x) d m 

, dx V-t^ — <i> (x, a) for m > , 

2 7 to! <9a m v ; 

dx /i(x)ef ] (a; \,\,T) + l [f dxdy h(x)4 £) (x y : A )h ! y ) 



(12a) 

(12b) 
(12c) 
(12d) 



Note that H^) describes the interaction of the scalar field with two flat plates and does not depend on the profile 



der in the profile. In we include counter terms of all orders in the profile h that depend 



h(x). H-^ is of m th or 

on the regularization parameter s but not on the quantum field <p. The constant c^'(a;A, A,T) of the one-point 
counter term depends on the plate separation a, temperature T and both coupling constants A, A. This finite counter 
term enforces the constraint of Eq.(|2]) at any temperature and separation when the interaction with the scalar <f> is 
turned on. It ensures that the parameter a represents the mean separation of the plates even when A and A do not 
vanish. [We shall later see that a different definition of the separation is to be preferred.] The coefficient function 
C2 (x — y; A) of the two-point counter term guarantees that the (measured) correlation (h(x)h(y)) at temperature 

T = remains Z?2(x — y) when the two plates are far apart and A > 0. c| (x — y; A) by construction does not depend 
on the separation a, temperature T or the coupling strength A of the distant plate. The (n > l)-point counter terms 
ensure that the corresponding 71-point correlation of the profile also remains unchanged at T = when the plates are 
far apart and the interaction with the scalar is switched on. These counter terms do not depend on a, T, orX and are 
oblivious to the existence of another plate, but diverge in the limit e — > 0+ in which the regularization is removed. 
The model requires an infinite number of counter terms and is not renormalizable in the sense of Dyson. However, 
it is renormalizable in a more modern sense (23j. Determination of the counterterms in fact does not diminish the 
predictive power of this model, since we assumed from the outset that correlation functions of the profile h at T = 
are all known (measured) when the other plate is far removed. No counter terms for the quantum field are required 
and the model unambiguously predicts finite effects due to interaction with another plate for any given profile h(x) 
[and its associated correlations] measured at T = 0. 

Thermal correlation functions of a free massless scalar field in equilibrium at temperature T in Matsubara's 
formalism [21] are generated by 



Zo[j;T] = exp 



-^ (0) + !B^ G nUn) 



(13) 



where J 7 ^ = — - ^ v is the Helmholtz free energy of a massless scalar field in a 3-dimensional Euclidean space of 
volume V ancQ 



tin\G° n \j n ) := / d 3 x / d 3 yj n (x)G° n (x-y)j n (y) 



The free thermal Greens-function, 



satisfies the differential equation, 



D -2-rrnT\x-y\ 

An\x — y\ 



(14) 



(15) 



(£ 2 - V 2 )G°(f - y) = 8{x-y) , with £ n = 2nnT. 



(16) 



4 The scalar product {. . . } defined in Eq.JTJ differs in its function space from (...) given in Eg. (I14D . We notationally differentiate between 
three-dimensional vectors u and two-dimensional vectors u. 
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The generating function of thermal Greens- functions at temperature T of the interacting model is[2l| 



Z[j, h; T, a] = exp 



exp 



Zo\j;T] 



-1 
Y 



^] + EE4 m, [^i 



(17) 



(18) 



Here Z"[j;T, a] generates the thermal Green's functions of the scalar field in the presence of two fiat parallel plates 
separated by a distance a, 



zHj;T,a] = exp 



_I Vff (0) f— 1 



^o[j;T] 



(19) 



exp 



exp 



) 2 +A ^ n (x,o) 2 - 



.I^II(T;a,A,A) + ^(i„|G||i n ) 



The free-energy J 7 " of a massless scalar field in the presence of two semi-transparent parallel plates was obtained 
in(l7| and is reproduced in App. [X] The thermal Green's function, g|, of a scalar thermal mode in the presence of 
two flat parallel plates satisfies the partial differential equation, 

{£ - V 2 + \8{z - a) + X5(z))Gl (x - y, z, z') = 5{z - z')S{x - y) , with £„ = 2imT. (20) 
Exploiting transverse translational symmetry, Gn is expressed by the dimensionally reduced Green's function as, 



ik(x-y) 



(21) 



with k 2 = £ 2 +k 2 = (27rnT) 2 +k 2 . Inserting Eq.([2l|) in Eq.® gives the ordinary second order differential equation 
satisfied by gW (z, z'\ k), 



(k 2 - -— + \5(z - a) + XS(z))g^ (z, z'; k) = 5(z - z') 
dz z 



(22) 



The solution to Eq. ([22l) for physical boundary conditions is well known [l7ll24| and reproduced in Eq. (|Bl[) of App. IBl 
Note that the reduced Green's function g" for the n th Matsubara mode is a function of K n only. An exponential 
cutoff, 



N(z,z';k) -> g ll(e) (z,z';K) = e EK g" (z, z'; k) 



(23) 



thus simultaneously regularizes loop integrals and summations. This regularization manifestly preserves transverse 
translational invariance of the corresponding regularized generating functional Z( e '\j,h;T,a]. Our renormalization 
condition that corrections to the correlation functions of the profile vanish at T = for a — > oo also is invariant under 
transverse translations. This implies that the regularized partition functional Z^[0,h;T,a] is analytic in the profile 
h with coefficient functions that are invariant under transverse translations. We therefore can use Eg. (1101) to evaluate 
this functional of the profile and express the free energy as, 



F{T) = J"H(T)-T lim lnZ (e) [0, 0; T, a], 



(24) 



where 



Z( £ )[0,0;T,a] = [0, T, a]Z h [c 



a=0 



With a Gaussian model describing the roughness correlations, the (regularized) generating function in Eq.([2"33 
is, 



Z^\j,a;T,a] = 

= exp £5>&>[&, £])] exp [iH^ 2 |«} + | E„(^|G! (£) |j n )] 



(25) 
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The functional [a, j; T, a] generates the Green's functions of two interacting scalar fields <fi and h that are supported 
on Euclidean spaces of different dimension. The base-manifold of the Euclidean quantum field 4> is S\ x IZ3, whereas 
the roughness field h is only supported on a two-dimensional plane in IZ3. From a contemporary point of view this is 
a miniature brane- world in which a matter field h(x) is confined to the two-dimensional universe of a plane embedded 
in a four-dimensional Euclidean bulk-space supporting <p. Because interactions occur only on the two-dimensional 
brane, the model is reducible to an equivalent 2 + 1-dimensional field theory on the 'boundary' universe of the plane. 
Roughness effects first manifest themselves in two-loop contributions to the free energy of this dimensionally reduced 
model. 



A. Feynman rules 

Eq. (|25|) defines a perturbative expansion and the associated Feynman rules for the scattering matrix on the plane. 
It will be advantageous to derive these rules in transverse momentum space. The presence of a second (flat) plate 
leads to non-local parts of propagators that are exponentially suppressed for momenta ko > 1. They describe the 
back-scattering off the distant (flat) plate and inform of its presence. 

1 . Propagators 

The model on the plane has four propagators. In transverse momentum space they are given by Eqs. (|B3a[) . (IB3bl) . 
()B3cP of appendix [Bl and by the Fourier transform c?(k) of Z?2- On the two-dimensional plane, 0„(x, a) and ^<^„(x, a) 
are independent and distinct modes. Introducing their Fourier components, 

*/>„(k) :=y"dxe jkx 0„(x,a) and ^ n (k) := j dxe lkx ^„(x, a) , (26) 



the four non-vanishing propagators of the surface model in (two-dimensional) Fourier space are, 

A^ : 



(f/> n (k)?/>„(- 


-k))H 


— 9ob\ K n) + 900 ( K n) 


1 

A + 2k„ 


(ip n (k)%j; n (- 


-k))U 


— <7[)1 \ K n) — 9iQ \ K n) 




(V>„(k)V> n (- 


■k))ll 


= 9ll On) + 9ll («n) 


~~ ~ T 


(h(k)h( 


-k)) 


= d^\ Ka ) = j dxD 2 (x)e ikx 





^ — 2ac t , a 



2A r 



(27a) 
(27b) 
(27c) 
(27d) 



with 



A„ := 1 - t n i n e- 2K " a , t n := \ , t n := \ , K n := v/(2^nT) 2 + k 2 . (28) 

In Eq.(|27j) we have decomposed the propagators into separation-dependent, (s)oft parts that are exponentially 
suppressed for a«> 1 and (/)ast parts that remain for infinite separation a — > 00. [In the following this distinction 
is dropped when irrelevant]. Note that (700 and <?oi vanish in the strong coupling (Dirichlet) limit A — > 00, whereas 
correlations of the normal derivative on the surface described by gn do not. There are no transitions between thermal 
modes in this model and the quantities t n , t n and A n defined in Eq.(|2"8]) are diagonal and functions of K n only. Lacking 
an experimental determination, the two-point correlation function for the profile is assumed to be given by a normal 
distribution in Eq. (|27dp . It is characterized by its variance a 2 and correlation length I. Although we will give results 
only for this particular form of the two-point correlation function, other correlation functions[13] that vanish faster 
than any power of the transverse momentum are equally admissible and do not qualitatively alter our considerations 
and conclusions. 

The dependence of the propagators on the exponential cutoff introduced in Eq. (|2"3"|) has been suppressed in Eq. (|2"7| 
but should be implicitly assumed. The cutoff length e can be neglected compared to the separation a in all expo- 
nentially suppressed (s)oft terms and loop integrals containing such a propagator remain finite in the limit e — > + . 
However, a-independent parts of internal loops containing only (/)ast propagators do diverge and a regularization is 
formally necessary to evaluate counter terms. 
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We collect the ip, ip propagators of Eg. ((27)) in the matrix, 



S(K)=[ 9 „°±\ 9 „°±\ ) . (29) 



9w{k) 9u («) 

whose negative inverse is the matrix of two-point vertices for vanishing profile h = 

A + 2k(1 + te~ 2Ka ) 2te~ 2Ka 



r^(/c):=-g («) = -( 2 f e - 2Ka -^(l-te- 2KO ) ' ^ 



-k \ / A \ - 2rea / k 1 

K- 1 j ~ I j ~ 216 [ 1 K" 1 



The corresponding generating functional of tree-level two-point vertices is, 



^n(-k) 
^n(-k) 



(31) 



(32) 



The dependence of r^ ^ on the coupling A is only in the -(^-component and is linear. A finite effective action implies 
vanishing ip but does not constrain ip in the strong coupling (Dirichlet) limit. Note that the quadratic form 1^°) [ip, ip] 
is an indefinite metric on the function space. 

The vertex function r(°)(K„) is diagonal in the Fourier-space of (k, £„)-modes and 

5 ln[- det r(°) ( K „)] = i ln(-^-) = \ ln(A„) + \ ln(l + + In 2 . (33) 

Comparing with Eqs. (|A1|) and (|A6|) shows that Eq. ([33l) essentially is the contribution to the free energy of a thermal 
mode in the presence of two parallel flat plates. Eg. (13"3l includes the contribution to the free energy due to the plate 
itself but not that of the other (distant) plate. This correspondence is further evidence that the (negative) effective 
action for the surface modes of a flat plate is indeed given by Eq. (|3"!?l) . 

2. Vertices 

Because the interaction in Eq. (TT2"j) is quadratic in the scalar <p and the profile h(x) does not depend on time, primitive 
vertices are diagonal in the Matsubara frequency and we need only specify their dependence on K n = \J f 2 + k 2 and 

The interaction H^,) in Eq. (|12c[) leads to transitions between the ip and '(/'-modes. It corresponds to the three-point 
vertex 7^ in fig. [T] 

7&W)=7£V, = (34) 

Expressions for primitive (m + 2)-point vertices j^ 171 ^ with m external roughness profiles are similarly obtained 
from Eq. (|12cp by noting that in Fourier space (<9/9a)" l 0„(k, a) may be replaced by k 2 ((9/9a)" l_2 0„(k, a) due 
to Eq. (|22p . Primitive vertices with an odd number of profiles lead to transitions between ip and ip modes and 
are given by, 

y§r 1 \n,*)=T$r l) V,K) = -X± ( 2n 2 t . (35) 

Vertices with an even number of profiles do not cause transitions between ip and ip fields and come in two kinds, 



7< 2 n W') = -A£ *? W a <»-*> and ^W) = -A£ U IM^^^ • ( 36 ) 

fc=0 ^ ' k=l ^ ' 



In diagrammatic form these vertices are shown in fig. [TJ 
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propagators: 



K„ 

— • 



■ ■+ «WW\^ 

¥„ goo V„ t g M ¥ ¥ En ¥„ h d h 



vertices: 




counter terms: 




-^(5")K 2k K' 2(n - k) 



X^AA t/l^Aj 



(e) 




V s 

- y£ d) - 

X £(2k n i) K2(k " 1)K ' 2< "" k) 



etc. 



FIG. 1. Propagators, vertices and counter terms of the 2 + 1 dimensional field theory on the planar surface. The 'roughness field' 
h corresponds to wavy- and the two dynamical surface fields to solid- and dashed- lines. Counter term vertices are depicted as 
crosses. Apart from ci, the theory only requires counter terms with an even number of external /i-legs. See the main text for 
details. 



Introducing the Fourier-transform h(k) = J e?xe lkx /i(x) of the profile, the primitive vertices may be collected to 
vertex functionals generating the interactions of the n th Matsubara mode with the profile, 



AVT(k,k'):=£ 



(2tt) 2 
; (2m)! 



2m+l 



6(k + k' + ^ kj) Joo n \ K n,K' n 



3=1 



2m + l 

5(k + k'+ ki)7S W+1) («nX) 
J'=l 



I/„ 10 (k,k'):=F n 01 (k',k) 



(37) 



Af n u (k,k') :=£ 



< (2m)! 



2m 



n 



dkj 



■ft(kj) 



2m 



<5(k + k' + ^k,) 7 ^ m) K,< 



Together with Eq. (f30l) the interactions of Eq. (|37| determine the vertex functional T[ip, ip; h] for any given profile h(x), 



^n(k') 



/ yOO ^01 \ 

where V n [/i] = " 10 t/ii 

\ "n n / 



(38) 



3. Counter terms 

The 2-point correlation function of the roughness profile of Eq. (127dp decays exponentially at large momenta and 
the vertex functional given by Eq. (|38[) is quadratic in the field <f>. One-particle-irreducible (1PI) vertex functions with 
only external ^-fields thus are finite if all 1PI vertices with only external roughness fields are. One therefore only 
requires counter terms for n-point vertex functions of the roughness profile. The 1-point counter term, cf is finite 
for e — > + and vanishes for a — > oo. This countcrtcrm is necessary for an unambiguous definition of the separation 
a. It ensures that Eq.Q holds at all temperatures, separations and couplings. The parameter a otherwise would not 
always represent the mean separation of the plates. Cjj is the only counter term that depends on the plate separation 
a, temperature T and both coupling constants A and A. To leading order in the loop expansion, the equation (h) = 
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FIG. 2. Feynman graphs for ci and C2 counter-terms to one loop, ci is finite but eliminates all tadpole contributions and 
guarantees that (h) = for any coupling, temperature and separation. Counter terms C2,ca,... are local and guarantee that 
corrections to prescribed roughness correlations vanish at T = in the limit a — > 0. 



is shown diagrammatically in fig. [5] and gives, 



(a; A, A) 



1-loop 



= -ta£ 



dk 



r#io( K n) 



~da~47 ^ 



/cd« ln(A) 



2ir\n\T 



-^/ (2) (r;A,A,a) 



(39) 
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where /' 2 '(T; A, A, a) is the Casimir pressure at finite temperature on two semitransparent plates due to a massless 
scalar field of Eq. (|A6[> . The last line in Eq.(l39| reproduces the Casimir force on Dirichlet plates at finite^] and at 

zero temperature [![. It is no coincidence that is the Casimir pressure since this counter term compensates for 
changes in the Casimir free energy due to (h) ^ 0. The correspondence is evidence for the correctness of the surface 

theory. Since it maintains (h) = 0, the counterterm cf~ ^ cancels all one-particle reducible contributions to the free 
energy. 

Counter terms with more than one external roughness field ensure that prescribed correlation functions of the profile 
remain unchanged at T — when the two plates are (infinitely) far apart. These counter terms by definition depend 
only on the coupling A and on the cutoff e and can be computed using the fast parts of propagators in Eq.j27j) that 

survive the a — > oo limit. Since g 01 = 0, counter terms with an odd number of external /i-fields vanish in the limit 

a — > oo. Apart from c\ the model requires only counter terms with an even number of external roughness profiles. 

To leading order in the loop expansion, (?; A) is obtained by evaluating the diagrams of fig [2] at T — in the limit 
a — » 0. For 1/A ^> e — > + one obtains, 



dk (k — n'e eK )e 



32tt 2 
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3A 3A 2 
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6e 



A 3 (l-31n2) 5gA 2 
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(40) 
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For Dirichlet boundary conditions one must consider the strong coupling limit < 1/A <C e — > + . The two-point 
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counter term in this case simplifies to, 



-2 (?;°o) = 2 



2tt 



dk 

(2^ 



32tt 2 



45 



g 

6e 2 



24 



(41) 



Together with the counter terms, the Feynman rules derived above define the loop expansion of this model. The total 
transverse momentum and thermal mode number are conserved at each vertex (assigning the time-independent /i-field 
the Matsubara frequency = 0). This is a 2 + 1-dimensional thermal field theory: the presence of another plate in 
a third spatial dimension manifests itself in the non-local dependence of propagators on the length scale " a" . From 
the point of view of the two-dimensional brane, this length scale could as well represent the Compton wave length of 
a massive particle. The model on the surface is holographic in the sense of(l8j. 



IV. THE DIRICHLET (STRONG COUPLING) LIMIT 



The vertices in Eqs. Q35[) and Q36p are all proportional to A. To leading order in the strong coupling expansion, 
the propagators <7oo>Soi,0io are of order A -1 and gn is of order A . The leading superficial order, N\, of a Feynman 
diagram in the strong coupling regime thus is given by, 



N\ = #700 + #7oi + #7io + #7n - #0oo - #0oi - #0io 



(42) 



where #JT denotes the number of Jf's the diagram is composed of. The quadratic model conserves the number of 
scalar surface fields and the ip and ip fields of propagators correspond to those of vertices. Vacuum diagrams thus 
satisfy the additional constraints, 



2#7oo 
2#7ii 



- #7oi 

- #7oi 
#7oi 



#7io 
#7io 
#7io 



2#<?oo 
#5oi = 



-#0oi 

- #001 
#010 ■ 



#010 
#010 



(43) 



Using Eq. (|43[) in Eq. (j42[) , the leading superficial order of a vacuum diagram is found to be given by the number of 
(?n propagators it contains, 



iV A (vac) = #3n = #7n + #goo - #7oo 



(44) 



If the strong coupling (Dirichlet) limit of the free energy is to exist, superficially divergent contributions with 
iV^(vac) > have to cancel. Such delicate cancellations generally arise due to underlying symmetries and are the 
consequence of associated Ward identities. A finite strong coupling limit for any profile in this sense is a non-trivial 
condition on the surface model defined by Eq. (I3"51) . That a Ward- like identity may ensure the existence of the strong 
coupling limit is suggested by the vertex functional for a flat plate given in Eg. (|3"2"1) . It evidently satisfies the 
identity, 







Sip n (k) d\ 



r<°) = o 



(45) 



Eq. (1451) can be interpreted as stating that for vanishing profile the normal derivative ip need not vanish when Dirichlet 
boundary conditions are enforced. The original interaction with the profile by the <5-function potential with the surface 
of Eq. (Ilip constrains the 0-field at strong coupling but not its normal derivative. The strong coupling limit otherwise 
would not correspond to Dirichlet boundary conditions. Even for non-vanishing profile, when ip and ip are coupled by 
V m [h] , the strong coupling limit must not require both surface fields to vanish. One therefore expects an /i-dependent 
linear combination of ip and ip to survive strong coupling and a generalization of Eq.(|32p to hold for the vertex 
functional T[ip, ip; h]. Writing the linear combination of thermal modes in terms of an /i-dependent functional 
the generalization of Eq. (l4"5l takes the form, 



<ty«(k) 



A n (k,k';h) 



(27T) 



<ty n (k') 



^■T[iP,iP;h}=0 



Inserting Eq. 



in Eq.pi)]) and varying tp(k) and ip(k) leads to the two functional relations, 
A n [h] ■ (1 - V™{h}) = V™[h] and A n \h] ■ V^[h] + V?[h] = . 



(46) 



(47) 
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FIG. 3. Cancellation of the leading order in A of contributions to ^-irreducible tpip- vertices. The solid lines represent goo- 
propagators which at strong coupling are A" 1 + 0(1). ^-irreducible ipip- vert ices thus are of leading superficial 0(A). For the 
vertices of Eqs. (|36[) and (|35|) . the leading superficial order cancels for any set of momenta and ^-irreducible iptf)- vertices in fact 
are of 0(1). 



A solution A n [h] to Eq. (|4"Tj) exists only if, 

V?[h] + V™[h](l - V™[h])- 1 V£ 1 [h] = , (48) 

for any profile h. We have explicitly verified Eg. (14"81 to sixth order in the profile h(k). Although we here do not 
provide a (combinatoric) proof of Eq. (|4"8")l to all orders, note that Eq. p5)l would determine A[h] and V 1 [h] for any 
choice of V 10 and V 00 . Requiring that solutions to the wave equation with Dirichlet boundary conditions are not 
trivial determines the interaction V^ 1 [/z,] in terms of y 10 [ft.] and y 00 [/i]. Eq. (j4"8|) could also be considered to generate 
7n-vertices that are consistent with proposed 701 and 700 vertices. It is an indication of the consistency of the 
model that vertices with up to six external roughness fields satisfy Eq. (|48|) . We will not require higher vertices in our 
calculations and may safely assume that Eq. (|4"5)) in fact holds to all orders. 

We still need to show that Eq.(|4"8]) is sufficient for a finite strong coupling limit of the effective action. In the 
following a connected Feynman diagram is ^-reducible if it becomes disjoint by removing a single <?n propagator and 
any number of d-propagators^. In this quadratic model, a vertex is ^-irreducible only if it contains no internal gn 
propagators. The analog of Eq.(|43p for a ip- irreducible ipift vertex diagram with two external ^-lines and no internal 
gn propagators implies that, 

2#7oo + #7oi + #7io = 2#.g o + #<?oi + #.9io 
2#7n + #7oi + #7io = 2 + #501 + #<?io ■ 
Its leading superficial order in A therefore is, 

iVA^-irred. ipip- vertex) = #700 + #701 + #710 + #7n - #.9oo - #.9oi - #.9io 

= #7n + #.9oo - #7oo = 1 ■ (49) 

Eq. (|48p on the other hand implies that the leading order in A of all contributions to an ■0-irreducible ^i/)-vertex in 
fact cancels. The superficial order in A in Eq. (|4"9")) does not account for this cancellation among contributions of the 
same superficial order and a V'-irreducible ^-vertex therefore is at most of order A . 

The superficial order in the coupling A of a vacuum diagram was found to be just #gn in Eq. (|4"4"]) . because this is 
precisely the number of ^-irreducible tp^- vertices the diagram contains. Since we have just seen that Eq. (l48p implies 
that a i/'-hreducible vertex in fact contributes at most is 0(1), the combined contribution to the free energy of 
all vacuum diagrams with a given number of <?n propagators also is at most of order 0(1). Eq. (|46[) thus ensures a 
finite free energy in the strong coupling (Dirichlet) limit. 



V. TWO-LOOP CONTRIBUTION TO THE FREE ENERGY: THE LEADING ROUGHNESS 

CORRECTION 

The two-loop vacuum diagrams of Fig. [4] give the leading roughness correction to the free-energy. The evaluation 
simplifies and is more transparent in the Dirichlet limit for both plates. The correction to the Casimir free energy per 



A diagram that can only be separated by cutting god 301 1 910 propagators and any number of h-lines is t/>-irreducible. A one-particle 
reducible diagram thus can be i/)-irreducible. 



13 




FIG. 4. Two-loop vacuum diagrams 



unit area of a massless scalar field for two parallel plates due to the roughness of one of them in this strong coupling 
limit is given by the rather simple expression, 



^' £ ' a ' T) = - T ^J (2,r)* (e^ - 1)(1 - e-*K) 



(50) 



where d(k) is the two-point correlation function of the roughness profile. There in addition is a finite correction to 
the free energy due to the roughness for an individual plate. It does not depend on the separation a and therefore 
does not lead to a modification of the force on the plate and will be ignored here. The correction to the interaction 
of Eq. ([50| depends on the exact form of d(k), but some conclusions about its general behavior can be drawn in the 
limit of large, £ 3> a, and of small, £ <C a, correlation length. For I ^> a, the support of the roughness correlation 
d(k — k'), such as the one of Eq. (|27d[) , is restricted to |k — k'|a <C 1. One thus may replace k' by k in the integrand 
without great error. With a 2 = (2ir)~ 2 J dkd(k) this gives the universal limit, 

Afg\a,l^a,T)~~ (JL\ /< 2 )(T; oo, oo, a) (51) 



-a/n 
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'} I-*.') LXJ 

~ Yd^ 2 2^ ^ Wna) 2 + (ni/T) 2 } 2 2T«<i' 240a 5 ' 

m=— oo n— 1 Lv 7 

which does not depend on the specific form of the correlation function d(k). As should be expected (25[, Eq. (|51"j) 
coincides with the roughness correction in PFA for profiles with large correlation length[l(| [HI- In the opposite 
limit of short correlation length, I -C a, or at large separations, the k-integral is exponentially restricted to the 
domain |k| < 1/a whereas the k'-integral in Eg. (1501) is finite only because roughness correlations are negligible for 
|k'| 3> l/£^> 1/a. The leading behavior of the roughness correction at separations a ^ £ thus is, 



A/^(i«o,A~cx>)~- 



J (2tt) 2 1 ; 4-7 (2tt): 



(e 2aK ™ - 1) 

, x ^/«(T;oo,oo,a) (52) 



= /W (T; oo, oo, af ff ) - /( 2 ) (T; oo, oo, a) + O , 
where f {2) {T; oo, oo, a) is the free energy of Ea. (IA6p for two flat parallel Dirichlet planes at separation a and 

fl ^ = a_ /(lr^ fcd(k) ~ a_ T\/! ' (53) 

The shift away from the mean of the profile is always of order a 2 /£, but the proportionality constant depends somewhat 
on the shape of the correlation function d(k) and is ~ 1-25 ... for the one of Eq. (|27dp only. Note that this 

displacement in the apparent surface of the profile is within the "thickness" of the profile for a < £. This mild 
condition generally is satisfied by naturally rough surfaces that arise from random dislocations of surface atoms and 
is a requirement for the validity of the loop expansion. However, it should be noted that surfaces with a > £ can be 
artificially created. In this case a loop expansion of the free energy in a 2 j£ 2 is not applicable ptl |27| . 
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FIG. 5. (color online) Relative roughness corrections to the Casimir energy and Casimir force in % due to a scalar satisfying 
Dirichlet boundary conditions on two plates, one of which is flat, the profile of the other is characterized by its variance 
a 2 — 49nm 2 and correlation length I. In two-loop approximation the correction is proportional to a 2 . Pairs of dashed and solid 
curves of the same color correspond to the same I — lOnm (violet), 15nm (blue),20nm (cyan),25nm (green) and £ = oo (from 
the outer pair of curves to the inner). Dashed curves represent the correction as a function of the mean separation o, whereas 
solid curves show it as a function of the effective separation a e ff = a — Zj- \f\- The (red) PFA correction for i — oo is the same 
in both cases. 



The perturbative roughness correction is well knownj^, 0, [j| [T^, [Hj|. A systematic loop expansion that includes 
temperature and roughness corrections simultaneously and yields a relatively simple closed expression as in Eq. (|50[) 
to our knowledge was not considered previously. Although it may not be worth the effort for the scalar field theory, 
the systematic inclusion of higher orders in the loop expansion in principle is straightforward in this field theoretic 
approach. 

Although the shift in Eq. (j53"|) generally is quite small and well within the profile's thickness, the effect on the 
roughness correction can be dramatic. As shown in Fig.[5l or as can be deduced by examining Eq. (|50[) . the perturbative 
roughness correction tends to increase with decreasing correlation length when the mean separation between the two 
plates is used as reference. The correction is quite large even for a ^> I and easily exceeds 20% at experimentally 
accessible separations for typical roughness profilesjsj. However, this effect to a great extent is eliminated by redefining 
the effective planar 'surface' of a rough plate. As shown in Fig. [5J the residual roughness correction decreases with 
decreasing correlation length I if the effective separation a® s of Eq. ([53|) is used for the separation. Thus, at least to 
leading order in the loop expansion, the main effect of roughness is to define the reference plane of the plate. This 
reference plane generally is not the mean of the profile. 

To determine the absolute separation of rough plates can be experimentally challenging. The previous considerations 
suggest that one could instead experimentally calibrate the (effective) separation of two plates so as to eliminate 
asymptotic 1/a 4 corrections to the Casimir interaction energy of fiat parallel plates (or asymptotic l/a 5 -corrections 
to the force). In terms of this more appropriate definition of the separation, the leading asymptotic correction to the 
force for large a £ is of order 1/a 6 only. Note that PFA-corrections, corresponding to infinite correlation length, 
are of this order and are not changed by this procedure. Such an intrinsic determination of the effective absolute 
separation a® s eliminates systematic errors due to electrostatic and other means of deducing the average separation 
of rough surfaces and facilitates a theoretical interpretation of the experiments. However, the suggested asymptotic 
calibration suffers from the fact that the Casimir force and therefore the signal-to-noise ratio decrease rapidly with 
increasing separation. A truly asymptotic determination is impractical and a compromise is required. Fig. [5] suggests 
that intermediate separations (lOOnm < a < 300nm) could be used to optimize this procedure in most experimental 
situations. In terms of the thus optimized separation, corrections to the Casimir force of two flat plates are much 
smaller and under better theoretical control. 

As argued in the introduction, an improved definition of the effective separation is necessary to avoid the conclusion 
that the reflection coefficient of a rough plate at long wavelengths (a ^> I) is larger than for a perfectly reflecting mirror 
and that unitary is violated. With a correct definition, the scattering matrix of a rough plate and the corresponding 
Casimir force ought to both decrease in magnitude compared to those of a flat (perfectly reflecting) Dirichlet plate. 
One furthermore expects the scattering matrix and Casimir force to decrease in magnitude with decreasing £ for 
a»f. Both physical requirements are met for a 3> £ by using the effective separation a^ s defined in Eq. (|53[) . The 
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FIG. 6. Localization of connected vertices in the I — ^ limit. Only (/)ast components of internal dynamical surface propagators 
contribute. External momenta are limited to |k| < l/a<l/f. 

corresponding force is always weaker than the PFA suggests. We now show that an appropiate definition of the 
absolute separation to a rough plate leads to a scattering matrix with physically acceptable properties in the limit 
£ > a. 



VI. THE LIMIT a > I: AN EFFECTIVE LOW-ENERGY FIELD THEORY. 

In the limit £ <C a the two-point correlation function /^(x) of the profile is localized to |x| < I <C a and we can 
use renormalization group techniques 28] to analyze the situation. In this limit we can approximately "integrate out" 
high momentum contributions and construct an effective theory of surface fields for wave numbers |k| < 1/a. In our 
model, the separation of momentum scales already occurs for tree- level surface field propagators. In Eq. (f2"Tj) they 
naturally decompose into (/)ast and (s)oft components. 

A local vertex v N ^ of the effective low-energy theory corresponds to the sum of all connected diagrams that contain 
only (/)ast internal propagators and have 2N (s)oft external ip- and 2N (s)oft external V'-fields and no external h- 
fields. Since the range of the roughness correlation Z?2(x — y) vanishes, these vertices are local in the limit £ — > and 
by construction do not depend on the presence of another plate at separation a and coupling A. Because g\^ = 0, 
the effective local vertices furthermore do not mix dynamical surface fields and conserve the number of ip- and ip- 
fields individually. They vanish unless the number of external ip and ip fields are both even. Closed loops of fast 
surface fields correspond to separation-independent corrections to roughness-correlations that are precisely canceled 
by the corresponding counter-terms at T — 0. At low temperatures it therefore suffices to consider connected vertex 
diagrams with only fast internal propagators and no closed internal loops of dynamical surface field propagators. The 
effective local vertices are finite since all loop momenta are restricted to \k\£ < 1 by the roughness-correlations, but 
some tend to diverge for I 0. To determine the degree of divergence with I, note that the number, L, of transverse 
momentum loops of a connected vertex diagram with 2(N + N) external surface fields is given by, 

L = #d+1-N-N (54) 

where #d is the number of roughness-propagators d(k) the diagram contains. The canonical mass-dimension [i^jv] 
of a N N- vertex in transverse momentum space is, 

[\~ 2N v Nfj ] = 2 - 37V - 3N . (55) 

In Eq. (|55l) a factor of A^ 1 provided by and propagators in vacuum diagrams was included for each external 
^-field. We distinguish two extreme limits: 

A. A rough Dirichlet plate: 1/A « f « a 

For 1/A <C £ the internal gj^ propagators (given in Eq. (|27aj) ) of an effective local vertex may be approximated by 
1/A. X£ 3> 1 includes the case of Dirichlet boundary conditions on the rough plate and we for simplicity consider only 
this extreme limit. The leading contribution to an NN- vertex is proportional to A^ due to the cancellations observed 
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in Sec. IIVI Note that internal d-propagators of a diagram are proportional to a 2 . For vanishing external momenta 
Eqs. ([5T]) and ([55)1 then imply that the sum of L-loop contributions to the effective local vertex behaves as, 

a- 2 ^-t^) JV+ " 1 © L " 1 ( 1 + °(^")) • (56) 

in the strong coupling limit. Upon summing the ^L-loop contributions, the local effective vertex at vanishing external 
momentum in the Dirichlet limit thus is of the form, 



where the dimensionless functions Q^^(s) are analytic at s = 0. We emphasize that the effective vertices reflect 
properties of the rough plate only. They do not depend on properties of the other parallel plate and we have the 
desired separation of scales. a 2 ji 2 < 0.1 for typical surfaces used in Casimir studics[8, 9] and low orders in the loop 
expansion should provide fairly accurate local vertices v^ti m * ne s t ron g coupling limit. 

Before proceeding to evaluate local effective vertices to leading order in the (hard) loop expansion, observe that the 
function in Eq. (|57[) depends on the two-point correlation function e?(k) and, in principle, also depend on higher 

correlation functions of the roughness profile. It therefore largely is a matter of perspective whether the vertices v N ft 
of Eq. ([57]) or the correlation functions D n of Eq. ^ are considered as the primary parameters that describe the plates 
roughness in the low energy effective field theory. Of course, not every set of local vertices corresponds to a 



physically realizable profile. A model of the correlations provides a basis for appropriate values and relations among 
the effective vertices of the low-energy description. Nevertheless, within a certain domain, the phenomenological 
parameters of the low-energy effective theory in effect are its local vertices. Assuming an analytic continuation of 
the functions Q^^(s) to s > 1 to exist, this effective low-energy description in the limit a ^ I can be extended to a 
region of the parameter space where the loop-expansion is no longer valid [26l [27j. 

Observe that the dependence of effective local vertices on external momenta of the (s)oft fields gives rise to con- 
tributions to the free energy that are suppressed by powers of ij a. To leading order in the loop expansion in a 2 / i 2 , 
Eq. (|57[) implies that effective local vertices with more than four (s)oft external fields can be ignored in the limit 
t/a — > 0. As in chiral perturbation theory[29] one arrives at an expansion in the canonical dimension of local vertices, 
vertices with more external fields becoming relevant at higher orders of the (soft) loop expansion only. Allowing for at 
most one hard internal loop, the present low-energy effective model (with the correlation function d(k) of Eq. (|27dp ) 
has the following local effective vertices, 
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Only the final expressions of Eq. ([58l) include the Kronecker symbols giving the conservation of Matsubara modes 
[indices with a bar designate ip modes] . Note that the tree-level (" one-roughness-exchange" ) and one- loop contributions 
to vn differ in their flow of mode indices. The individual terms in the intermediate expressions correspond to diagrams 
of different topology that contribute to the effective vertex. 

The free energy of the effective low-energy theory with the (s)oft propagators of Eq.([2T]) and local vertices of Eq. ([55|) 
describes separation-dependent corrections due to the profile of a Dirichlet plate at separations from a smooth 

parallel plate. The free energy of the low-energy effective theory in powers of a /a may be obtained in the (s)oft loop 
expansion. The effective 2-point vertices i>oi and t>io play a crucial role in this expansion. They correct the low-energy 
behavior of propagators and thus affect all higher orders of the (s)low loop expansion as well. 
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FIG. 7. One loop contributions to the free energy in the effective model for a rough Dirichlet plate. 



The ratio v® : — A 2 /4 precisely compensates for the ratios of soft propagators g[f : g[ S Q : g^tf — —A/ (2t) in the 
Dirichlet (t — > to = 1) limit. The roughness and separation-dependent correction, A/ 1 - 2 ', to the free energy per unit 
area of the effective low energy model thus is obtained by evaluating the 1-loop diagrams of Fig. [7] with an effective 
interaction —2p D = and propagator g[f , 
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= f E / 7^2 [ ln (! ~ Cgh(«n)^e- 2 ^«?«) - ln(A„)] (59) 

n 

The one-loop free energy depends on the mean plate separation a and on a length p D = —vqi ~ V2na 2 /£ that 
characterizes the profile. The effective separation of the two plates, = a — p D /2, in one-hard-loop approximation 
coincides with the one found perturbatively in Eg. (I5U1) . We in addition obtain the reduced scattering matrix t^ugh 
for low-energy scattering off a rough Dirichlet plate, 

CW = d + AK K . (60) 

^rough( K ) ^ S positive and never exceeds unity. It satisfies all the requirements of a reduced scattering matrix and is 
consistent with phenomenology for scattering off a rough plate in that only short wavelengths with np D 3> 1 are 
strongly affected. ^ughC' 1 ) < 1 i s due to diffuse scattering of part of the incident wave with (transverse) wave vector 
k. The intensity of the outgoing wave with (transverse) wave-vector k is thereby reduced. Diffuse scattering by 
a rough surface is more effective at shorter wavelengths and negligible for wavelengths that are long compared to 
p D a 2 /I. It perhaps is worth remarking that the scatterin g m atrix found in this approximation does not depend 
on the separation a of the two plates (as the GTGT- formula|16j indeed requires). However, our approximations in 
deriving the low energy effective theory can only be justified for a»f 



B. A rough semi-transparent plate: a>f <1/A 

This limit includes that of weak coupling. We proceed similarly as for the Dirichlet case but I now is the smallest 

correlation length. The leading behavior of a local vertex thus is determined by its degree of divergence as the 'cutoff' 

( f \ ( f \ 

£ on hard loop momenta is removed. In the limit of large transverse momenta we have that gQ ~ l/k,g\i ~ —k/2. 
Neither depends on A and a local vertex in this case is proportional to X Nv , where Ny is the total number of primitive 
vertices it is composed of. Eqs. (|55|) and (|54|) then imply that for vanishing external momenta, 

\2jV . /„2\# L ~ 1 
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The meaning of Eq. (f6Tj) becomes more transparent upon noting that for any connected diagram Ny — 2N > N—N+l. 
At any given order in the loop expansion, the largest contribution to an effective local vertex in the present limit arises 
from diagrams with the minimal number Ny = N + N + 1 of internal vertices. Local n-point vertices with external 
V^-ficlds thus are suppressed by powers of A£ <C 1 compared to vertices that describe the scattering of slow -0-quanta 

(s) 

only. In the present limit, the low energy effective model is described by a scalar surface field with propagator <?q 
and local vno vertices only. To first order in the loop expansion one again obtains the vertices of Eqs. (|58bj) and (|58c[) . 



poo i \ 2 2 

v 10 = -X 2 £ 2 a 2 kdk-e- k2£2 / 2 > - ^-V^5 nn , (62a) 

J 2 A£«i U 

v 2Q = 2^V / kdk2 ^ e~ fc2£2 > ^±^ 5nn , 5mm , , (62b) 

Jo V 2 / A£ «! 2 

(62c) 

but the local effective interactions v$\, v\\ and 1*02 now are negligible. To one (s)oft loop, the roughness and separation- 
dependent correction A/( 2 ) to the free energy per unit area in the limit a 3> I -C 1/A is, 

(ik x ^ 1 ( ^"nXn^n^ 71 



2 ^ J (2tt) 

with 



(2tt) 2 2k V A„ 
| E j (^2 N A » - P^ni„e~ 2K " Q ) - ln(A„)] 
| E / (^2 [ ln (! - (! + PKn)a„e- 2 ^) - In(A n )] 
r E / 7^2 N 1 - WK)^- 2 ™) - ln(A„)] , (63) 



irough('«) = (1 + np)t(n)e Kp < t(n) and a c ff = a — p/2 . (64) 

One reproduces the Dirichlet-plate result of Eq. (|60| by simply setting t = t D = 1 and p = p D in Eq. (IMl) , but for the 
same profile the parameter p is only half that found in the Dirichlct limit, 

The displacement of the equivalent surface of the rough plate from the mean of the profile by p evidently also depends 
on the transparency of the plate. Considering the effective shift p(X) as a phenomenological parameter of the rough 
plate, the main effect due to roughness in the limit a 3> I is to define the position of the planar scattering surface 
and simultaneously modify the scattering matrix of a flat plate as in Eq. (|64p . It is interesting that the modification 
of the scattering matrix and the shift from the mean of the profile are not independent. The two extreme limits we 
considered provide a range for the parameter p(X) in terms of the variance and correlation length of a profile described 
by Eq.(j22d|, 

< P(\) < ■ (66) 

The upper bound of Eq. (|66[) corresponds to a rough surface with Dirichlet boundary conditions and the lower to 
weak coupling. Note that the effective scattering plane does not coincide with the mean of the profile even for weak 
coupling A ~ [although the scattering matrix is arbitrary small]. 



VII. DISCUSSION 



We have developed a field theoretical approach to the Casimir free energy of a massless scalar field in the presence 
of parallel rough and smooth semi-transparent plates. Changes in the free energy due to interaction of the scalar with 
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FIG. 8. (color online) Relative roughness corrections to the Casimir energy in % due to a scalar satisfying Dirichlet boundary 
conditions on two plates, one of which is flat. The profile of the other is characterized by its variance a 2 = 49nm 2 and 
correlation length I. The separation is between equivalent planes representing the plates (see the text and Eq. (153p for its 
relation to the mean separation.) The leading two-loop approximation for different correlation lengths t is given by solid 
curves that correspond to those of fig. [5] Dashed curves represent the correction in the effective low energy theory derived 
in the limit a 3> I. Pairs of dashed and solid curves of the same color correspond to the same correlation length I = 
lOnm (violet), 15nm (blue),20nm (cyan),25nm (green) and I = oo (red). The leading two-loop approximation interpolates 
between the low-energy model for large separations a ^> I and the PFA result (solid red) for small separations a < I. Note that 
typical roughness corrections are much smaller than the PFA suggests. 



the rough surface were found to be described by an effective 2 + 1-dimensional held theory on the equivalent plane 
involving two dynamical surface helds, and tp as well as the static prohle h. The model on this planar boundary 
of the original space is holographic in that the existence of another dimension and of a second parallel plate at a 
separation a are encoded in non-local propagators. The theory in this sense is a low-dimensional analog of branc 
models in string theory [l8j. 

Two-loop contributions to the free energy of this model give the leading roughness correction. This leading correction 
for a massless scalar field is qualitatively similar to that obtained by perturbative analysis for electromagnetic helds [5, 
S 0, HH , but its field theoretic origin allows for a consistent inclusion of finite temperature effects and for a more 
transparent interpretation. For a scalar field in the strong coupling (Dirichlet) limit, the leading loop correction was 
obtained in closed form in Eq. (|50]) and is shown in fig. [5j As for the electrodynamic corrections considered in [9l. ll2T] 
the PFA result [l(J [H| is reproduced for a<< and the Casimir force appears to strengthen with decreasing £/a. We 
argued in the introduction that this apparently violates unitarity when a>l 

The problem can be traced to an inappropriate choice of the equivalent planar surface for a rough plate. This 
plane does not coincide with the mean of the profile but is displaced a distance p oc a 2 j t from it. With the correct 
definition of the effective surface, roughness corrections are much smaller and the Casimir force weakens with increasing 
roughness a 2 /£. Roughness strengthens the Casimir force only for a /£ < 0.5 and only for a < At. In this regime the 
unitarity argument based on transverse translational symmetry does not hold. In terms of the effective separation, 
the PFA to the roughness correction is approached from below with increasing correlation length. As pointed out at 
the end of Sect. [V] it should be possible to intrinsically calibrate the absolute separation using experimental results 
and take advantage of the much smaller corrections. 

We finally derived an effective low energy field theory that describes the limit a £ with a single length parameter 
p ~ a 2 /£ characterizing the roughness of a plate. The correction in this limit indeed is described by an effective 
scattering matrix trough for a plane displaced a distance p/2 from the mean of the profile as given in Eq. (164p . As 
illustrated by fig. [5J roughness in the effective low energy theory weakens the force at all separations and the effective 
scattering matrix is always less than for a flat plate of the same material, approaching the scattering matrix of the flat 
plate for wavelengths k < 1/p. It is also evident from fig.[5]that the corrected 2-loop estimate interpolates between the 
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low energy effective model and the PFA, approaching the former for small and the latter for large correlation length 
£. For realistic correlation lengths and variances of the profile, the roughness correction at most is a few percent at 
a = lOOnm for a scalar field satisfying Dirichlet boundary conditions. It is even less for semi-transparent materials. 

At separations where roughness corrections are important, temperature corrections are small and vice versa. The 
numerical results for small separations shown in figs. [5] & [8] therefore do not show temperature corrections. However, 
the formalism and most of the equations we derived include them. In the experimentally relevant electromagnetic 
case this could be of interest in the transition region 500nm-2/i where both corrections are small but comparable. 
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Appendix A: Free Energy of a Massless Scalar Field for Two Flat Parallel Semitransparent Plates 

1. An isolated flat semi-transparent plate 

Although this contribution to the free energy does not depend on the separation a of two flat plates, it is finite and 
does depend on the temperature. We compute it for the sake of completeness. 

Using Matsubara's formalism one[24| readily finds that the irreducible contribution to the Helmholtz free energy 
per unit area, /W, of a massless scalar field due to a semi-transparent flat plate of area A described by the potential 
interaction V(z) — X5(z) is given by, 



/W(T,A) 



T 



E 



n— — oc 



dk 



(27T) 



2 Mi 



A), 

2hi n 



(Al) 



where T is the temperature and k 2 — (2imT) 2 + k 2 . Poisson's resummation formula allows one to rewrite Eq. (jAl[) in 
the form, 



n=-oo J ~°° J v ' 

00 rp pOO \ 

= E^7 o ^sMWT)ln(l + -) 
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where the divergent, but temperature- independent, n — summand has been dropped by requiring that the free 
energy vanish at T — 0. This amounts to ignoring the divergent change in zero-point energy due to insertion of 
a semi-transparent plate. In deriving the second expression of Eq. (IA2[) we introduced spherical coordinates with 



C, 2 + k 2 and performed the angular integrations. The final expression in Eq. (|A2p is in fact finite. We may 
perform the summation and reduce the expression for the free energy per unit area of a flat plate to a single integral, 
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The asymptotic behavior of f^ 1 ' is readily found, 

/«(T«A)~ — 
/ (1) (A«T) 
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(A4) 
(A5) 
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For Dirichlet boundary conditions (A — > oo), the asymptotic expression in Eq. (|A4l) holds at any tempera- 
ture. Eq. (|A5p is accurate to leading order in A for a weakly interacting plate. Note that the free energy of a 
single semi-transparent plate is positive and increases monotonic with temperature for any value of A. The corre- 
sponding contribution to the entropy therefore decreases with increasing temperature. However, this ignores the bulk 
contribution to the entropy which generally overwhelms this reduction. Including the bulk contribution, the total 

entropy due to insertion of a Dirichlet plate is negative only for 1/T > 13 5 V/A ~ 2V/A. It is negative only when 
the boundary of the container (on average) is within a thermal wavelength of the plate. Ignoring the finite size of 
the container in obtaining the entropy change due to the plate is no longer warranted in this situation. Although we 
here do not quantify the correction, it very likely is perfectly consistent that the entropy change due to insertion of 
a single plate is negative and decreases as the temperature increases. The negative contribution to the entropy can 
be qualitatively understood by the fact that the energy difference for excited cavity states increases upon insertion of 
the plate and the occupation numbers for excited states therefore decrease. 



2. Irreducible contribution to the free energy of a scalar due to two flat parallel semi-transparent plates 

We again use Matsubara's formalism and proceed as for a single plate. The irreducible contribution to the free 
energy per unit area, f^ 2 \ due to two semi-transparent parallel plates at separation a is given by, 

/( 2 >(T;A,A,a)^ £] / ^ ln(A( K „)) = ^ £ / K d«ln(A(«)) , (A6) 



TL — — O0 TI— — 0O 



where k 2 = (2nnT) 2 + k 2 as before and A(k) is given by Eq. (|B2l) . Contrary to the irreducible contribution from a 
single plate, is finite for any separation a > 0. We again use Poisson's resummation formula to express the free 
energy in dual variables, 

-^rH^P^h™ (A7> 

cLkkN( ) In 1 



2tt J y 2nT' V (A + 2k)(A + 2k) 

Here N(x) is the staircase function ([x] denoting the largest integer less than x), 



N(x) := 1/2 + [x]=x + - arctan(cot(7rx)) . (A8) 
At low temperatures behaves as, 

1 f°° 7T 2 T 4 

/( 2 )(2^7a« 1;A,A) ~ — / dn n 2 ln(A) + Aa——— , (A9) 
4tt 2 J 90 

where the effective separation a = a + j + ~. The first term is just the Casimir energy of two semi-transparent plates 
[13] ■ Note that the T 4 behavior of the second term is the same as that of the bulk contribution to the free energy. In 
the Dirichlet limit A, A ~ oo it simply subtracts the contribution to the free energy from the volume between the two 
plates. This again is qualitatively caused by the increased energy difference to excited states between the plates. The 
second term in Eq. (|A9p is not correct in the weak coupling limit when 2irT S> A, A. In the range A, A <C 2ttT <l/a 
we have that 

/ (2) (A,A< 2^T< I /a) - (A10) 

Xl ^ ■ n n / ^ln(T/A)-Aln(T/A) , i WMa \ 19 2 
~ 32^ V + 2 " Ta { A - A + L27036 J ~ U {2nTa) + • ' ■ 

Note that for weak coupling the entropy apparently diverges like ln(T) for small T. However, there is no violation of 
Nernst's theorem in this case, because Eq. (|A10[) only holds for 2irT ^> A, A. For lower temperatures Eq. (|A9p is valid 
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and the entropy vanishes proportional to T 3 . The first term of Eq. (|A10[) reproduces the leading term of the Casimir 
energy for two weakly interacting parallel plates [13, [HI- 

The total free energy ]F\\ , of a massless scalar field in the presence of two parallel flat plates is the sum of the bulk 
contribution, the irreducible one-body contributions of the individual plates in Eq. (|A3p and the irreducible two-body 
contribution of Eq. (|A7[) , 

J"" (T; A, A, a) = + AfW (T; A, A, a) , (All) 

with 

J" (T; A, A, a) = /W(T, A) + J«(T, A) + f^(T; A, A, a) . (A12) 

We have chosen the normalization of the free energy so that J-"" vanishes at T = for widely separated plates, thus 
absorbing a divergent, but temperature and separation independent, factor in the normalization of the generating 
function. 



Appendix B: Thermal Green's Function of a Scalar in the Presence of Two Parallel Semitransparent Plates 



In Matsubara's formalism [2 1| thermal Green's functions of a mode at temperature T are given by evaluating 
Euclidean Green's functions at the corresponding Matsubara frequency £ n = 2nnT. We thus can draw on the literature 
for the Euclidean Green's function of a massless scalar in the presence of parallel semitransparent plates [lil Il7l. |24| . 
The physical solution to Eq. (|2"2"|) is, 
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Of particular interest to us is the correlation function in momentum space at z = z' = a and its derivatives (</4(x, o) 
^</>„(x,a),<(x,a) = ^r0„(x,a), etc.), 
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where the expression is to be evaluated at the n-th Matsubara frequency (k — > K n — y/ (2irnT) 2 + k 2 ). The 
correlations in Eq. (|B3l) are found by taking normal derivatives of Eq. (|Bip and using that lim s _>o sign(s) = 0, 
lim^o sign 2 (s) = 1 and lim s ^o<5(s) = lim s _ ! .o sig n '( s ) = 0. Eq. (|B3el) expresses the fact that Eq. (j2"2")l relates 
correlations on the surface of the rough plate to ones with two fewer normal derivatives of <p. Increasing the num- 
ber of normal derivatives by two amounts to multiplying the Fourier-space correlation function by k 2 . The three 
correlation functions of Eqs. (|B3al) . (|B3b|) and (|B3cp thus generate all correlations with a higher number of normal 
derivatives such as Eq. (|B3d[) . This allows us to obtain Feynman rules for vertices with an arbitrary number of /i-fields 
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